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Abstract 

A state sum construction on closed manifolds a la Kuperberg can be used to construct the partition 
functions of 3D lattice gauge theories based on involutory Hopf algebras, A, of which the group algebras, 
CG, are a particular case. Transfer matrices can be obtained by carrying out this construction on a manifold 
with boundary. Various Hamiltonians of physical interest can be obtained from these transfer matrices 
by playing around with the parameters the transfer matrix is a function of. The 2D quantum double 
Hamiltonians of Kitaev can be obtained from such transfer matrices for specific values of these parameters. 
A initial study of such models has been carried out in [T]. In this paper we study other regions of this 
parameter space to obtain some new and known models. The new model comprise of Hamiltonians which 
“partially” confine the excitations of the quantum double Hamiltonians which are usually deconfined. The 
state sum construction allows for parameters depending on the position in obtaining the transfer matrices 
and thus it is natural to expect disordered Hamiltonians from them. Thus one set of known models consist 
of the disordered quantum double Hamiltonians. Finally we obtain quantum double Hamiltonians perturbed 
by magnetic fields which have been considered earlier in the literature to study the stability of topological 
order to perturbations. 


1 Introduction 

Topological Order - Topologically ordered systems have gained wide attention in recent years due to some 
of its consequences in topological quantum computation and emergence of new phases of matter among many 
others [5] . Among the different types of systems exhibiting topological order the ones with long-ranged entangled 
(LRE) ground states are the ones which are thought to be most useful for quantum computation. The earliest 
proposals of such systems are the quantum double Hamiltonians of Kitaev [311]. The toric code is the simplest 
example of a 2D lattice systems which contains anyons as low energy excitations and have degenerate LRE 
states as ground states, this model consists of spin variables living on the links of a square lattice. These 
were further generalized by the Levin-Wen models (3 or the string-net models which described more general 
anyonic excitations by directly taking a unitary fusion category as inputs. These models are also quantum 
double models based on weak Hopf algebras as noted in [B] and can thus be constructed via the algorithm of 
Kitaev. Topological codes have also been considered on manifolds with boundary [ZlIH]. Several other models 
inspired by the usefulness of the toric code as a stabilizer code have been constructed of which the topological 
color codes [3iini[n] are an example which have also been experimentally implemented |12| . 

State Sum Constructions and Statistical Mechanical Models - Both the Levin-Wen model and Toric code 
models can be thought of as Hamiltonian realizations of topological field theories (TQFTs) and then they can 
be formulated in terms of topological invariants. State sum constructions of TQFT’s |13j have been employed in 
realizing statistical mechanical models in the past [HI [TS] [16] . Such methods have also been used to construct 
the Levin-Wen models [S] using the Turaev-Viro invariants m and chain-mail link invariants [HIH]. Kitaev’s 
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toric code has also been related to Turaev-Viro codes [^D]. The Levin-Wen model corresponds to a topological 
invariant called Barrett-Westbury invariant m and the toric code corresponds to a special case of the Kuperberg 
invariant [I]. This has been especially noted in [22]. We showed this explicitly in [I] where we embedded the 
2D quantum double models based on an involutory Hopf algebra, A, in an enlarged parameter space (defined 
later), that of the 3D generalized lattice gauge theories based on these algebras A. In the special case where 
the algebra A is taken to be the group algebra CG of a group G, the generalized gauge theory can reproduce 
the lattice gauge theories familiar to physicists US). The toric code occurs when we choose G = Z 2 [I|. 

The Quantum Double Model - The quantum double model of a discrete group gQ is defined on a bidimensional 
lattice over a compact manifold E of genus g. The degrees of freedom live on the links of the lattice and they 
are vectors \g)i G TLi, where g € G and I represents a link of the lattice. The total Hilbert space TL is then the 
tensor product of all Hi (for all 1), in other words H = Hi®H 2 ® ■■ ■ ® Hn, where N is the total number of 
links. A basis vector |^') of H is then of the form |'I') = |( 7 i)i (g) 132)2 ® • • • \gN)N, with gi G G. The dynamics 
of such a model is governed by a Hamiltonian made up of a sum of commuting operators, acting locally 

on the plaquettes, p, and vertices, v of the lattice and this is given by 

= -JpY,Bp-JvJ2Av , ( 1 ) 

p V 

where Bp is the plaquette operator, the vertex operator and Jp^y positive numbers. These operators are 
both projectors and also commute with each other for all vertex and plaquettes making energy levels discrete. 

Consider {|a) : a = 1,2, • ■ • , fc} a complete set of eigenvectors of with |a) = Ea\a). If the system is 
put into a bath with temperature T one can obtain thermodynamics properties of such a system by its partition 
function 

Z^«=^(a|e-'5^|a)=tr(e-'^^), (2) 

a—1 

with P = l/ksT, being ks the Boltzmann constant. The matrix is called the transfer matrix. Since H is 

given by equation (TJ the partition function can also be written as 

Z^^{P,Jp,Jy) . (3) 

\ p V / 

In the special case where = 0 this partition function can be rewritten as 

Z^‘3(/3,Jp,0)=tr (n ^hJpBp 

) = En g/3JpS(p) ^ 

\ P / conf. P 

where the sum runs over all the configurations and S{p) = -1-1 if the holonomy of the plaquette p is flat and 
B{p) = “1 otherwise. The function S{p) is invariant under gauge transformation, and then Z^^{P, Jp, 0) is the 
partition function of a lattice gauge model. For Jy ^ 0 it can also be shown that Z^^ is a partition function 
of some lattice gauge model [Tj. 

Partition Functions of 3D Lattice Gauge Theories - We mentioned earlier that the quantum double Hamilto¬ 
nian is related to the partition function of a lattice gauge model, as once we know the Hamiltonian the partition 
function is well defined. Now we can ask the question, if a partition function Z is given, is it possible to obtain 
a Hamiltonian H such that the equation [2] is satisfied? The answer for this question is no, however if Z is 
restricted to be gauge invariant it can be done. The reason is that there is a way of building gauge invariant 
partition functions out of a 3D topological invariant called the Kuperberg invariant |24| . which is based on 
involutory Hopf algebras. Moreover this construction allows us to obtain partition functions that are more 
general than the lattice gauge ones, but they are still gauge invariant. 

This partition function is parametrized by four non-physical parameters, namely zs,Zt G CC? and G 

CG*, in other words Z will be a function of the form Z{zs,Zt,^sAt)- The choice of such parameters leads to 
specific models. It is important to note that these parameters are not free such as for example /3 in ordinary 

^The quantum double model is defined for more general inputs like involutory Hopf algebras and more generally weak Hopf 
algebras. We will only be concerned with the case of group algebras which are a particular case of involutory Hopf algebras. 
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gauge theories. These parameters are fixed once the model is fixed. In particular, as it can be seen in [53], if 
the parameters are taken to be H 


Zs = Zt = 


2® 


V+i 


2 ® 


(5) 


and = e, where e denotes the counit of the algebra CG. The partition function obtained by such a 

choice coincides with that of a 3D pure lattice gauge theory with G = 2 - Different choices of these parameters 
lead to different partition functions which may not be related with gauge theories at all, however if we restrict 
the parameters zs and zt to be elements of the centre of CG, whatever the partition function is, it will be gauge 
invariant |23j . The action in this partition function may not be of physical interest but they will all be gauge 
invariant, and for that reason we will say that Z{zs, zt,^St^t) is the partition function of a generalized gauge 
theory. From now on we consider G to be a discrete group, and the 3D dimensional lattice a finite cubic lattice 
as a triangulation of a manifold of the form S x where S is a 2D compact manifold of genus g and is 
the one dimensional sphere. 

Kuperberg’s Construction of Transfer Matrices - The way to build these generalized partition functions 
Z{zs, zt,^S,^t) is by associating tensors, made up of the structure constants of the algebra, to the faces and 
links of the lattice leading to a very complicated tensor network that, fortunately, can be realized as the trace 
of a matrix U{zs,zt,^s,^t) by 


Z{zs,zt,^sAt) =tl'{U{zs,ZT,^S,^T)) , ( 6 ) 

where, by analogy, we can think of it as being the exponential of a Hamiltonian H{zs, zt,^sAt)j i- e. 

Uizs, ZT, ^S, ^t) = , (7) 

where the /3 constant can be suppressed without loss of generality. The matrix e~^ can be thought of as an 
operator acting on a 2D lattice (over the 2D manifold E), but for that we have to make distinction between 
the timelike and spacelike directions on the original 3D latticcH. This procedure is shown in detail in [1] , where 
a Hilbert space H is associated with the 2D lattice as 77 = "Hi 0 H 2 0 • • ■ Hat, with Hi ^ CG being the local 
Hilbert space associated with the link I of the 2D lattice and N is the total number of links of the 2D lattice. 
The degrees of freedom are then elements of CG living on the links, which is equivalent to saying that the 
degrees of freedom are group elements, since there is a one-to-one correspondence between group elements and 
elements of the basis of CG. 

Models from the transfer matrix- We then obtain the Hamiltonian by taking the logarithm of U{zs, zt, ^s, ^t)- 
However the Hilbert space H is very huge, which makes U{zs, zt,^s, Ct) a very huge matrix and difficult to take 
its logarithm. But as we have done before [I] the matrix U{zs, zt, ^s, Ct) can be decomposed into a product of 
local operators acting on H given by 

U{zs,zt,^s,^t) = 

p I V 

where p, v and I denotes plaquettes, vertices and links respectively with T/ and L/ operators acting on the 
degrees of freedom located on the links I and Bp and Ay the plaquette and vertex operator which we previously 
encountered in the quantum double Hamiltonian given in Eq. |TJ These vertex and plaquette operators satisfy 
the quantum double algebra Moreover the plaquette and vertex operators commute with each other 

for all choices of plaquettes and vertices. The link operators do not commute in general and thus switching 
on the parameters corresponding to them namely, zt^^s will complicate the procedure of taking logarithms of 
the transfer matrix. Therefore we can not obtain exactly solvable Hamiltonians for an arbitrary choice of the 
parameters 2 : 5 , zt, ^nd we can only do it for those of which the local operators commute with each 

^Here is the basis for the group algebra CZ 2 . 

^The terms timelike and spacelike are just used to distinguish the spacelike directions in 2D from the third direction which we 
call timelike. We still work in the Euclidean metric. 

^The notion of quantum doubles arises in the theory of Hopf algebras where the quantum double construction is used to generate 
a quasitriangular Hopf algebra from a given Hopf algebra. A quasitriangular Hopf algebra is governed by a i? matrix which satisfies 
the quantum Yang-Baxter equation (QBYE). This can also be taken to be a way to generate solutions for the (QBYE) which 
coincide with the irreducible representations of the braid group in two dimensions. In physical terms these are anyons which are 
also the irreducible representations of the quantum double algebra and hence the usefulness of the quantum double Hamiltonians in 
obtaining anyons in the spectrum of the theory. The reader is referred to the book on Hopf algebras where these ideas are discussed 

m- 
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other. One such Hamiltonian is the quantum double Hamiltonian of Eq. [T] where only the parameters zs and 

corresponding to the plaquette and vertex operators are used. 

Other models using more parameters in the transfer matrix - Some examples include quasi-topological phases 
which result from a condensation of the excitations [27] of the quantum double phase of Kitaev. This leads 
to increased ground state degeneracy for the condensed phases. Examples of these phases were studied in [T]. 
These phases including the quantum double phases of Kitaev were obtained when we considered the parameters 
Zs and f,T in the transfer matrix. In m we showed that we could obtain the quantum double phases of Kitaev 
by writing down models which included zt or ^ 5 . 

Identifying topologically ordered phases - The models that exhibit topological order in 2D can be understood 
by their quasi-particle content called anyons. The data that determines the phase are the ground state degen¬ 
eracy, their statistics and fusion parameters [52]. If two models have different ground state degeneracies, or 
different fusion rules or statistics, they are not in the same topological phase. 

In this paper we take this program further by considering more parameters in the transfer matrix which were 
not included in [11 [28]. We consider three types of models here. Two of them do not include zt and fs while 
the third includes them. The first two sets of models comprise of the disordered quantum double Hamiltonians 
of Kitaev and a new Hamiltonian which leads to “partial” confinement of the excitations of the quantum double 
phase of Kitaev. We use the term partial to emphasize the fact that the models are such that the excitations can 
be moved a few steps with an energy cost after which they become deconfined like in the usual quantum double 
models. According to the terms added we can confine the excitations for any number of steps that we wish to. 
We will also call these models n-step confined models in the text to follow. The models which include the other 
two parameters, zt and fs are the quantum double Hamiltonians perturbed by magnetic fields. These can be 
thought of as local perturbations to the exactly solvable Hamiltonians of Kitaev. Due to the usefulness of these 
models to realize fault tolerant quantum computation, it is necessary to study the stability of the topological 
order to local perturbations [201 Ell Eg. These models have already been considered in the literature and we 
write them down here just for the sake of completion and to drive home the point that they are well within the 
parameter space of the three dimensional lattice gauge theories. Our focus is on exactly solvable Hamiltonians 
like the original toric code Hamiltonian and so phase transitions are out of the scope of this paper as we will 
then necessarily have to move through perturbed toric code Hamiltonians which are outside the exactly solvable 
regime. 

The contents of the paper are organized as follows. Section 2 gives a brief review of the construction of 
the transfer matrix of the generalized lattice gauge theories. The section also includes an introduction to the 
mathematical preliminaries that go into the construction of the partition function and the transfer matrices. The 
algebra of operators, which include the quantum double relations between the vertex and plaquette operators, 
are written down. The models obtained from this transfer matrix are described in section 3. An outlook is 
presented in section 4. 


2 Partition Function and Transfer Matrix of Generalized Lattice 
Gauge Theories 


The partition function of lattice gauge theory is a well known example of a classical partition function built 
out of local weights associated to plaquettes of an oriented 3D lattice, where the gauge degrees of freedom are 
elements of a gauge group G living on the edges of the lattice. A configuration is a choice of an element g € G 
for each link of the lattice. For G = 2 the gauge degrees of freedom are spin variables ±1 living on the links. 

The action of this model is defined by S' = ^ where the sum runs over the plaquettes of 

the lattice and Up is the holonomy of a plaquette p. The partition function which describes the model is given 

by 

Z = ^ Y[M{Up) ,with M{Up) = exp{-/3/2 (tiiUp) + . (9) 

conf. conf. P 

In above equation M{Up) is the local weight for the model. Due to the invariance of the local action under cyclic 
permutation, the local weight is also invariant under this cyclic permutation, which makes it a class function, 
M : G —>■ C (equivalently, M{g) = M{hgh~^)). This construction can be generalized by choosing M{g) as being 
any class function M : G ^ C. Moreover, we can associate local weights, A(/), to the edges I of the lattice, 


4 


( 10 ) 


such that the partition function is now given by 

conf. P I 

We can reproduce the usual lattice gauge theories by making appropriate choices for M{p) and A(Z) and also 
generate partition functions that are still gauge invariant but do not represent a physical model. The partition 
function in equation ra is called the partition function of a generalized lattice gauge theory. Starting from 
this partition function we can obtain a transfer matrix whose logarithm gives us Hamiltonian operators, and 
thus dynamical quantum models defined over the Hilbert space % defined before. These quantum models are 
parametrized by functions of the parameters of the generalized lattice gauge theories. In [T] it was shown that 
the quantum double Hamiltonians of Kitaev, of which the toric code is a special case, can be obtained from this 
approach. In other words it was shown how to embed such models in the parameter space of these generalized 
lattice gauge theories. 

In [T] it was shown that this partition function can be build out of the structure constants of an involutory 
Hopf algebra A and a 3-manifold of the form E x 5'^, where E is some compact 2-manifold and is the 
1-dimensional sphere. We did not consider all possible deformations of the generalized lattice gauge theory 
partition function in [1], working only with a specific kind of deformation (one parameter deformation) of the 
gauge theory partition function. Now we will allow other deformations by letting the parameters be any element 
of the center of the algebra A and it’s dual algebra, A*. We only work with group algebras CG of a discrete 
group G here. Nevertheless the methods presented here hold for any involutory Hopf algebra. We will go 
through the mathematical preliminaries beginning with the definition of the group algebra CG. 

2.1 The Group Algebra CG 

The group algebra CG of a discrete finite group G is generated by the basis elements {(jig : g G G} indexed 
by the group elements and let {vl® : g € G} be the dual basis such that = 5(g, h). In this basis the 

multiplication and co-multiplication are defined by 


4^a-4^h ■— 4^ab Pblab — Si^db^C) 

:= (5(a,6)V'“ ^ = S{a,c)6{b,c) 


and we can easily see that the unit and the co-unit of the algebra are 


e 

e 



where e is the identity element of the group. 
Finally the antipode map is defined by 


S {4>a) = 4>a~^ 


Ca = S(a,e) 
e“ = 1 


= S (ab, e). 


It is not difficult to see that the group algebra structure constants satisfy all the axioms of Hopf algebras [551 
[T] . An important thing about the group algebra is the fact that it is an involutory Hopf algebra, which means 
that = 1 . 

As the simplest example lets consider the group algebra of 2 - The group 2 is defined by 2 = {il}, and 
the product is the usual multiplication. The group algebra C 2 has {4>+\,4)-i\ as basis, and the product is 
given by 


(p+l4>+l — ; 

4>-i4>+i = 4'+i4>-i = 4>-i ■ 

The dual basis is and the coproduct defined by 

\]/+ivl/+i = v[/+i ^ 

= xjz+lxlf-l = 0 . 
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The antipode in this case is trivial S = 


2.2 Constructing a Partition Function with an Involutory Hopf Algebra 


Consider an oriented 3D cubic lattice as a triangulation of a closed 3D manifold of the form E x and let 
G be a discrete finite group. The degrees of freedom are elements of the gauge group G located on the links 
of this lattice. A configuration is a choice of one group element for each link of the lattice. The generalized 
partition function is then built by associating tensors for each face and link of the lattice which will play the 
role of local weights. These tensors contract with each other resulting in a scalar called generalized partition 
function Z(z, ^). 

For a plaquette, whose boundary links carry the group elements a, b, c and d, we associate a tensor Mated 
defined by the algebra structure as shown in figure |l(a)| For a link with four plaquettes glued to it’s edges, 
labelled by x, y, z and t, we associate a tensor defined by the coalgebra structure as shown in figure 1 1 (b) I 

It is very important to have an algebra and coalgebra structure so that the associated tensors can contract with 
each other on the closed manifold. These tensors are parametrized by the elements of the center of the group 


d 


a 



c 


Mabcd 



^xyzt 


(a) The tensor Mated associ¬ 
ated to a plaquette of the lat¬ 
tice. 


(b) The tensor asso¬ 

ciated to a link of the lat¬ 
tice. 


Figure 1: Local weights associated to the plaquettes and links of the lattice. 


algebra CG and it’s dual, z and ^ 

Mabcd{z) = tr (z (j)a4>b4>c4>d) , (11) 

and 

A“^^*(z*) = co-tr . ( 12 ) 

where tr((/)g) = \G\5{g,e) is the trace in the regular representation and co-tr('I'S) = 1 , V 5 . The partition 
function is obtained by contracting the indices of the tensors associated to the plaquettes and links. However 
we need to take care of the orientation of the lattice while performing this contraction. If the plaquette and 
link orientation matches the contraction is made directly, otherwise it is done through the antipode tensor S'^. 
At the end the partition function will be of the form 

Z{Z,0= E \[Mated{z)\{/\^y^\0\{Si , (13) 

indices P l !■' 


where the sum runs over all the contracted indices of all tensors and I' runs over the links with mismatching 
orientations. 

As an example take G = 2 and choose 2 :® = | (e^^+i + and = e. Then the tensors will be of 

the form 


Mated — {zy(f>ated') — 2^^ ((^^^+1 A ^ ^4^—l')4^ated') 

= I {e^tr{(j)ated) + e~^tT{(j)-abcd)) 

— Q/3{a.ted) . 

/^xyzt ^ co-tr = co-tr 

= S{x,y)S{x,z)6{x,t) . 


(14) 

(15) 


Note that abed is the holonomy Up of the plaquete p. In this case, due to Eq. [131 the sum over indices in the 
partition function in Eq. [TS] will reduce to a sum of indices of the tensor Mated (for all plaquettes), in other 
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words, it will become a sum over configurations and the partition function will take the form 

indices P l conf. P 

that is exactly the partition function of the lattice gauge theory defined in Eq. [HI 

We can now extract a transfer matrix U{z,^) out of this partition function, but first we have to make a 
distinction between timelike and spacelike directions of the 3D lattice, as shown in figure [2] for a small piece of 
the lattice. Since there is now a distinction between the timelike and spacelike parts there is nothing forcing the 



Figure 2: Time and spacelike direction. 


parameters of the timelike and spacelike plaquette weights to be the same. So for a more general description we 
have all spacelike plaquette weights parametrized by the element zs El while the timelike ones are parametrized 
by zt- In a similar manner the spacelike and timelike link weights are parametrized by £md respectively. 
Therefore the partition function is now a function of the form Z{zs, Zt, ^t) El- 

2.3 The Transfer Matrix 


From the partition function we have just defined we can get a transfer matrix U such that its trace is equal 
to the partition function. This transfer matrix may depend on the same parameters as the partition function, 
namely U = U{zs, The operator U acts on the links of the 2-dimensional lattice where the quantum 

states lives. A local Hilbert space Hi associated to each link I (with basis {|g)i : g G G}). The Hilbert space 
of the full system is given hy H = Hi ® H 2 ® ® Hn and a vector of this space is a linear combination of 
vectors of the form 

l5i)i ® 152)2 (D • • • 0 \gN)N ■ (17) 

The procedure to get such a transfer matrix was shown in [1] using a diagrammatic notation to manage the 
tensors that builds the partition function and also the transfer matrix. We are now considering a more general 
parametrization than the one considered in [T]. The way to get the transfer matrix is similar to the procedure 
shown in [T], resulting in the same operators as in [I] but now in a bigger parameter space. We just write down 
the results in what follows. The most general transfer matrix including all the parameters is given by 

U{zs,zt,^s,^t) = ]^-Sp(2s) ]^Gz(zt,^s) A„(^t) , (18) 

p l V 

where Bp(zs) is an operator which acts on the links at edge of the plaquette p, A„(^j-) an operator which acts 
on the links sharing the vertex v and Ci{zt,^s) is an operator which acts on a single link. The operators Bp{zs) 
and Ay{^T) are called the plaquette and vertex operators, as before. All the parameters in U{zs, Ztt^Sj^t) are 
central elements of A and A* . It can be seen in [33] that the central elements of CG are written in terms of the 
conjugacy classes [G] of G as 

2 ; = E , with zc = (t>g , (19) 

c ge[C] 

^The parameters zs need not be the same for all the spacelike plaquettes. This shows that this construction need not obey 
translational invariance. This fact will be exploited to obtain disordered Hamiltonians in the next section. The same fact holds for 
the other parameters as well. 

®The partition function also depends on the group G and the lattice ZZ, so it is actually a function of the form 
Z(G, C, zsi ZTiiSi^T)- For brevity we denote it just as Z(zs, 
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while the central elements of A* we write in terms of the irreducible representations i? of G as 

^ = E = E , (20) 

R geG 

where xr{9) is th® trace of g in the representation R. Hence the parameters of the transfer matrix can be 
written as 


= E 

C 

= E 

R 


and Zt = E 
C 

and = E 
R 


The plaquette and vertex operators are linear functions in its parameters. That is 


c 


BS 


Avi^r) = E ^R Ay (C'^) = E o-rAv ■ 


R 


where the operators B^ and A^ act on the links as shown in the picture in figure |3] and are given by 


( 21 ) 

( 22 ) 


(23) 

(24) 




V 5 


1 

k p > 

3 



2 




Figure 3: A plaquette and a vertex of the lattice. 


= |G| E <5(ar^a2a3a4\G)rfi0r“=®r3“^0r“S (25) 

{a.ltl 

Av = '^XR{9)R'i{g~^)®Ri{g~^)®Lb{g)®LQ{g) (26) 

g&G 

where S{g, G) = 1 if g G [G] and 6{g, G) = 0 if g ^ [G] and the operators Li{(j)g), Ri{(j)g) and T;)'?®) are operators 
which act on a single link defined as 


Lli4'g)\h)l = \gh)i , 

Rli(t^g)\h)l = \hg)i , 

Ti{^g)\h)i = S{g,h)\h)i. 


(27) 


These operators are linear on its parameters, in other words, L{z) = z^L[(j)g) (this property also holds for 

R{z) and T{z*). Sometimes we use the short notation L® = L{(l>g), i?® = R{4>g) and T® = r('I'®). The link 
operator Ci{ztAs) can also be written in terms of the Li and T/ operator as we shall see next. 

The main difference between this model in the one we have considered in [T] is the link operator Ci{zt,^s)- 
Unlike in [T] this operator is no longer proportional to identity, now it takes the form 


Ci{zT,is) = \G\ms)Li{zT) 


|G| 


E««^'(E) 

V ^ EE 


/ 


c ' ' 


ip 


(28) 
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where 


L? = R? = E = E 

se[G] se[G] 

Ti^ = T.XRi9)Tl> 


g£G 


Thus the final expression for the transfer matrix is 


u{zs, ZT, cs, ct )=\Gr n E n 


P \ c 


J2^RTn E^^^ 


c 


i L \ H 


n E 


« \ R 


(29) 

(30) 


R 

V > 


(31) 


where the parameters on the left hand side are related to the coefficients on the right hand side by the equations 
(Ell) and 1^. 

2.4 Algebra of the operators 


In order to find the algebra that the operators Bp, Ay and C/ satisfy one should first look at the algebra of 
the operators Li, Ri and T;. It is not difficult to see that the following relations holds 


Lfif = Lf i?f = d{g, h)Tf 

= RhL9 j^gj^h ^ rj^haj^g j^grj.h ^ j^hgj^g 


(32) 


Now using Ea. (l5^ and the orthogonality relations on the characters xr{9) [33] we can show the algebra of the 
operators which build the transfer matrix, namely Bp , Ay, Lp and Tj^. 

The set of operators {Bp} is a complete basis of orthogonal projectors which generate the plaquette opera¬ 
tors, which means 

B^B^' = \G\6{C, G')B^ , and IGp^ ^ = 1 • (33) 

G 

Same way the set of {A^} is a complete basis of orthogonal projectors which generate the vertex operator, in 
other words 

A^A^' =5{R,R')A^ , and ^4^ = 1. (34) 

R 

The plaquette and vertex operator still commute for any choice of the parameters zs and so we can write 


[Bp,Ay\=0 ^ [i?p(2s), A(^t)] = 0 , V2:s,^t- 


(35) 


The set of operators {Lf} and {T^-} are also complete sets of orthogonal projectors {Lf Lf = 6{G,C')Lf and 
= S{R, R')Tp'), however the operator does not commute with the link operator but it does commute 
with the plaquette operator, in the same way as the operator Lf does not commute with the plaquette operator 
but it does commute with the vertex operator. Thus we can write 


[Tf^,Bf] = Q, [r«,A«]^0, [Lf,Bf]^0, [Lf,Af]=Q. 


(36) 


Therefore, the link operator Gi{zt,^s) commutes with the plaquette and vertex operators only for some choices 
of the parameters zt and ^s, which means the quantum model obtained from the transfer matrix containing 
this link operator is in general not solvable. 


3 Examples of Models from the Transfer Matrix 

The fully parametrized transfer matrix obtained in the previous section helps us construct a number of other 
interesting models. We have seen that the quantum double Hamiltonians are only one special class of models 
in this parameter space. Here we will look at what other possibilities exist in the extended parameter space. 
The first set of examples consist of disordered quantum double Hamiltonians. These are quantum double 
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Hamiltonians which do not have translational invariance. This is due to the appearance of coefficients, for the 
vertex and plaquette terms, which are not constant but depend on the vertex v and plaquette p. These models 
continue to be in the quantum double phase. The trace of the transfer matrix for these Hamiltonians also help 
us obtain their partition functions. 

We can produce solvable models even in the presence of the parameters zt and ^ 5 . We discussed a class 
of such models in |28j where we exhibited exactly solvable models which continued to remain in the quantum 
double phase described by modified vertex operators or plaquette operators. Here we show another class of 
models that can be obtained from the transfer matrices of lattice gauge theories that continue to remain exactly 
solvable but the excitations are “partially” confined with respect to those of the original quantum double model. 
By this we mean that the deconfined quantum double excitations are now confined up to a few steps due to the 
addition of an extra term in the Hamiltonian. The number of steps for which they are confined can be controlled 
by adding the appropriate term in the Hamiltonian. We will call these models n-step confined models. These 
models comprise our second set of examples. They are examples where the self-duality of the quantum double 
Hamiltonians is broken. We also discuss the ground states, it’s degeneracy apart from the excited states of the 
model. 

Finally we write down models that are obtained by using the remaining two parameters zt and along 
with zs and ^t- The transfer matrix is now significantly modified as new single qudit operators, acting on 
individual links, appear. They do not commute with the vertex and plaquette operators in general. However 
for certain special values of parameters they commute with products of vertex and plaquette operators as we 
shall see when we consider these models later in this section. The models obtained at these values are quantum 
double Hamiltonians perturbed by generalized “magnetic” fields. They have the interpretation of magnetic 
fields in the case the input algebra A is £(^ 2 ). For the remaining values of the parameters we obtain more 
complicated terms in the Hamiltonian which we will briefly touch upon. These models are not exactly solvable 
and are outside the phase described by the quantum double Hamiltonians, at least for sufficiently large values 
of the parameters. 


3.1 Disordered Quantum Double Hamiltonians (QDH) 

The transfer matrices used to obtain these models only use zs and ^t- The other two parameters are set to 
Zt = T] and ^ for all the timelike plaquettes and spacelike links. To obtain the disordered QDH models we 
associate a different central element of the algebra and it’s dual to every spacelike plaquette and timelike link 
respectively. This leads to the following transfer matrix 

(37) 

p V 


where zs^p and ^t,v are the parameters for the plaquette p and vertex v respectively. The plaquette and vertex 
operators commute with each other in this transfer matrix and each of the operators is a sum of projectors. 
Thus it is easy to take the logarithm of these matrices to obtain the disordered QDHs. 

Let us look at these Hamiltonians in the case when A — €(^ 2 )- Denote the basis elements of €(^ 2 ) by 
with = (j)i, and the basis elements of the dual €(^ 2 )* by with the product = 

The central elements of C{Z 2 ) and €(^ 2 )* can be written as 2 ; = ai(j)i+a-i(j)-i and ^ = /3itp^+/3-itp~^ 
respectively. By assigning each spacelike plaquette and timelike link with zs,p and ^t,v respectively we obtain 
the transfer matrix as 

U {£{Z2),zs,pAt,v) =W_Bp{zs,p) (38) 

p V 

where the plaquette operators are given by 

Bp{zs,p) = oti^pBl -h a-TpB~^ (39) 


with 


B. 


±1 


and the vertex operators are given by 


1 ± al 


2 


(40) 


Av{^T,v) 




Pl,v 


- P-l,v 

2 


A. 


-1 


(41) 
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with 


(42) 


A4. . 


The action is shown in figure (O. 

The disordered QDH can now be written as 


H = 




+ 1,44 ^ + X! (inai.P-Sp + Ina-i^pB^ i) . 


(43) 


As it can be seen this Hamiltonian breaks translational invariance and is made up of a sum of commuting 
projectors. When translational invariance is restored we recover the familiar toric code Hamiltonian. 

The ground state for this Hamiltonian is easily obtained by projecting on to the smallest of and 
for every vertex v and ai^p and Q;-i,p for every plaquette p. The ground state degeneracy is the same as the 
usual toric code and the winding operators, 


Xct,cixci,ci) 


n 




(44) 


^ci.c2.(ci.c2) - n 

feeCi,C2,(Ci.C2) 

where the non-contractible loops Ci, ( 72 , Cj*, (7| are defined on the direct and dual lattice respectively as shown 
in the figured! commute with the Hamiltonian. 


^2* C2 



(a) Non-contractible loops on the direct 
and on the dual lattice. 



(b) Non-contractible loops on 
a torus. 


Figure 4: Non-contractible loops. 

The excitations correspond to the other value of the coefficients for each vertex and plaquette. The string 
operators (ribbon operators in the case of the non-Abelian groups [27]) creating the excitations are the same 
as in the QDH case. For the specific case of £(^ 2 ) we have the string operators, creating charge or vertex 
excitations at the end points of the string 7 along the direct lattice, as 

= (46) 

and those of the fluxes or plaquette excitations at the end points of the string 7 * along the dual lattice as 

Pt = n (47) 

fce7* 

These are shown in figure dl As is well known these excitations are deconfined by which we mean that there 
is no cost in energy for moving them around by stretching the string creating them. Moreover the fusion rules 
and braiding statistics are the same as in the translationally invariant toric code. Thus we conclude that the 
disordered QDH given in Eq. continues to remain in the toric code phase. 
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Figure 5: Deconfined flux excitations. 


3.2 Quantum Double Hamiltonian with n-Step Confined Excitations 

We first consider the case with n = 2. To obtain this we write down an example of exactly solvable Hamiltonian 
made up of the QDH vertex and plaquette operators along with new terms made of these operators which have 
2-step confined low energy excitations unlike the QDH case where all the excitations are completely deconfined. 
The transfer matrix used to obtain these models contain only the zs and parameters. 

The transfer matrix can be written as 


U {A,Zs,iT) =WBp{zs) 

p V 


(48) 


In the case of C(Z 2 ) we can write down the Hamiltonian with the 2-step confined charges and fluxes as follows 


H = 

+ 

+ 


(oiH), -I- a-iA^ -I- -|- P-iB^ 

V p 


-I- a'2A\.A^^ + a'^A^^A],. + a'^A^^A^^'^ 

<*i> 

E + P',b;}b;}) 




(49) 


where < ij > and < i*j* > are nearest neighbor vertices in the direct and dual lattices respectively. All the 
terms in this Hamiltonian commute with each and other and are sums of projectors. The ground states are 
given by the usual toric code Hamiltonian. The degeneracy does not change as the winding operators in the 
toric code case, given by Ea. (j44p and Eg. (1451) . continue to commute with this Hamiltonian and thus help create 
the new states from a given ground state. In particular on a torus the degeneracy is four. 

The interesting feature of this model occurs when we look at the excitations. As in the toric code case the 
string operators creating charge and flux excitations are given by Eq. (H51) and Eq. (HTl) respectively. However in 
this model when we create a charge or flux excitation we also excite the direct or the dual link given by the 
Al.Al. or the B^,B^. 0 term respectively. This creates link excitations along the string where the operator 
given by Ea. (H51) or Ea. dTf)) acts. However the creation of these link excitations occurs only for two steps after 
which they are deconfined as in the toric code case. Thus we say that these excitations as 2-step confined 
excitations. These partially confined excitations are shown in figure (© and ([7]). 

It is easy to see that the dyonic excitations are also confined in a similar manner. Thus we have a model 
based on lattice gauge theory which is exactly solvable, has ground state degeneracy and has excitations which 
are confined up to two steps or can be thought of as being partially confined. 

There is a natural way to increase the number of steps for which these particles are confined. This is achieved 
by coupling more number of vertex and plaquette operators. For example to obtain three-step confinement we 
add terms of the form A),. where f,j, fc are nearest neighbors on the direct lattice. The corresponding 

plaquette terms are BhtBhtBht where k* are nearest neighbor vertices on the dual lattice. The energy of 

Pi Pj Pk 

these excitations are more when compared to the original QDH as we violate more terms in the Hamiltonian 
to obtain them. This argument can easily be extended to any number of steps. The corresponding figures of 
the n-step confined excitations will have a larger shaded region where they are confined when compared to the 
n = 2 case. 

^ We assume that the coefficients of these two terms in the Hamiltonian are the smallest and hence violating these will lead to 
excitations. There is no loss of generality in making this assumption. 


12 













E 


E + AE 


E + AE 


E + AE 


(a) The pair of charges 
with energy E. 


(b) Moving by one 
step increases the en¬ 
ergy by AE. 


(c) Moving one more 
step outside the 
shaded region does 
not cost any energy. 


(d) Once outside 
the shaded region 
the charges are 
deconfined. 


Figure 6: The 2-step confined charge excitations. 



(a) The pair of fluxes 
with energy E. 



(b) Moving the flux by 
one step increases the 
energy by AE. 



(c) Once outside the 
shaded region they are 
deconfined. 


Figure 7: The 2-step confined flux excitations. 


This model can be constructed on any triangulation of a two dimensional lattice keeping the confinement 
properties unchanged. This is easy to see as the additional terms of pairing neighboring vertex and plaquette 
operators is independent of triangulation. The model can also be easily extended to all group algebras and more 
generally to all involutory Hopf algebras without any obstacle. 

In the case of C{Zn) we can also add the parameters zt and to include the single qudit terms in the 
transfer matrix and keep the above properties of confinement. We will illustrate why this is so for €(^ 2 ). The 
argument extends easily for other n. We have the following relation 

Al^A^EAl^Al = (50) 

where Ti = —and vi, V 2 are the end points of the link 1. This reduces the model to the previous Hamiltonian 
showing confinement. From this identity we see that we can include the parameters in th® transfer matrix 
without any effect to the properties of the Hamiltonian considered earlier. 

Also note that the braiding and fusion rules do not change in the n = 2 case. The only difference from the 
usual toric code is that the energy of the quasiparticle excitations is now higher and there is an energy cost to 
move them up to one step in the lattice. Since the topological data is the same we conclude that this model is 
in the same phase as the toric code. By increasing n we could enter a different phase, namely a confined phase, 
in the thermodynamic limit. 

3.3 Perturbed QDHs 

These are constructed out of transfer matrices which include the parameters Zt or ^s- Their inclusion introduces 
single qudit operators Li[zt) and Ti{^s) on the link I respectively. These terms do not commute in general 
with the vertex and plaquette operators thereby making the process of taking their logarithms and hence 
obtaining the Hamiltonian difficult. However for certain parameters we can still take the logarithm to obtain a 
Hamiltonian. We will illustrate this in the case oi A = €(^ 2 )- 
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Let us include the parameter zt into the transfer matrix. This brings in the operator 


+ (51) 

when Zt = xi(j)i + X-i4>-i. When Li{zt) = erf we can have a product of two plaquette operators adjacent to 
the link I commute with Li{zt)- We can then take the logarithms to obtain the following Hamiltonian 

H {ln{ai)Al + ln{a-i)A-^) + Y + in • (52) 

V p ^ 

This model resembles adding magnetic perturbations to the QDH for the simplest case of €(^ 2 ). The dual 
version of this involves using the parameter in th® transfer matrix to obtain the magnetic field operator af 
on the links. 


4 Outlook 

A systematic procedure was used to obtain the models with interesting properties which included studying all 
possible ways of taking the logarithms of the transfer matrices of generalized lattice gauge theories. In [T] this 
led to quasi-topological phases with increased ground state degeneracy occurring due to condensed excitations 
of the QDH. In this paper another way of taking the logarithm led to partially confined excitations. We went 
further to include other parameters in the transfer matrix which took us away from the topologically ordered 
phases. 

The state sum procedure used in [1] and further explored here can be thought of as a method to construct the 
quantum double of a given input algebra. Using this principle we can use other inputs to obtain the quantum 
doubles of more general objects leading to more interesting models. One such input are groupoid algebras 
which are examples of quantum groupoids. Such considerations lead to confined excitations in pure lattice 
gauge theories m- In [Q the quantum doubles of weak Hopf algebras were constructed. These reproduced 
the Levin-Wen models. Our considerations will embed the Levin-Wen models in the parameter space of lattice 
gauge theories based on these weak Hopf algebras. The weak Hopf algebras used in [B] were the ones constructed 
from a unitary fusion category by Kitaev and Kong in [52]. In particular by using weak Hopf algebras of [22] as 
inputs in the construction of [T] we will be able to obtain the operators creating excitations in the Levin-Wen 
models and finally we can confine and condense these quasi particles by using the methods in this paper. 

The state sum construction can be used to construct the transfer matrices of lattice theories with gauge 
and matter helds by adding matter degrees of freedom, acted upon by the gauge fields, on the vertices of the 
triangulated lattice. This construction was shown to produce exactly solvable quantum models in one and two 
dimensions in |35j . The methods of this paper when applied to the transfer matrix constructed in |35) is bound 
to give many new interesting models. 
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